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ABSTRACT 

The coagulation of dust aggregates occurs in various astrophysical environments. Each one is characterized by different conditions 
that influence the growth, e.g., relative velocities, composition, and size of the smallest constituents (monomers). Here we study the 
microphysics of collisions of dust aggregates in a four-dimensional parameter space. The parameters are the collision energy, the 
initial compactness of agglomerates, the mass ratio of collision partners, and the impact parameter. For this purpose we employ a 
state of the art molecular dynamics type of model that has been extensively and successfully tested against laboratory experiments. It 
simulates the motion of individual monomers interacting dynamically via van der Waals surface forces. The structure of aggregates 
is quantified by the filling factor that provides information about the internal structure, the packing density of monomers, and the 
projected surface area of aggregates. 

Our results show the importance of the impact parameter that causes formation of elongated particles, due to tensile forces acting in 
off'set collisions. In head-on impacts, aggregates are compacted at lower energies. A sufficiently high energy causes restructuring to 
reach maximum compaction. If more energy is provided, pancake-like structures are formed. We find that the outcome of collisions 
can be represented in a simple way. A highly pronounced large fragment component coexists with a power-law distribution of small 
fragments. The structural parameter of these small fragments is described very well by a simple relation, largely independent of the 
initial compactness, impact energy or impact parameter. The simulations show that erosion by collisions with high mass-ratio can be 
significant. The ejected mass can be several orders of magnitude higher than the impactor mass. This contrasts with collisions of equal 
mass aggregates, where the same impact energy can lead to perfect sticking. 

These findings are summarized in the form of a simple collision recipe. The recipe specifies the outcome of a collision, averaged over 
the impact parameter. It is provided in tabular form for a range of physical parameters such as impact energy and pre-collision filling 
factor. The dependence on the mass ratio of impactor and target is taken into account by providing both a local and a global branch of 
the recipe. 
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1. Introduction 

.1^ Dust plays an important role in many astrophysical environ- 
^ ments, and it is also the building material of planets. Small 
' grains can stic k to each other thanks to van der Waals forces 
(I Johnson et al.lll97ll: iDerjaguin et al.lll975l) . Depending on the 
growth mechanism, the aggregates can acquire different struc- 
tures. In low velocity collisions, which are dominant in the 
regime where particle velocities are dominated by Brownian 
motion, the narrow size distribution primarily leads to col- 
lisions between particl es of a similar size, in turn leading 
to fractal aggregates (iKempf et al.l 119991: iBlum et al.l 120001: 
iKjause&Bliimil 120041: iPaszun & Dominikll2006h . Fractals pro- 
duced by Br ownian growth have a typical fractal d imension 
of Df = 1.5 (iKrause & BlumI 120041: iPaszun & Domln ik 2006), 
and higher gas densities shorten the mean free path of parti- 
cles and result in even more fluffy structures. In the limiting 
case of very high densities, aggregates formed this way may 
be very elon gated with the fractal dimension approaching unity 
(IPaszun & D ominik 2006). 

Velocities induced by turbulent gas motions and by radial 
drift in protoplanetary disks can be much higher. Furthermore, 
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the dependence of velocities on particle size can become re- 
versed, with larger particles moving faster than smaller particles. 
For a discussion of the relative velocities in protoplanetary disks , 
we refer the reader to a review article bv lBeckwith et al.l (l2000h . 
Relative motion that emphasizes collisions between particles of 
different sizes leads to the production of aggregates with very 
different structures, because then the particles also grow by col- 
lecting smaller project iles. In this case, the results are porous , 
non-fractal aggregates (iBafl & Wittenlll984l:IPaszun et al.ll2009h . 

Low-impact velocities generally result in sticking. However, 
the growth of dust usually causes an increase in collision veloc- 
ities when the particle start to decouple from the surrounding 
gas. When the collision energy becomes higher than the energy 
needed to roll monomers over each other (later referred to as 
rolling energy £"^11), restructuring begins. Very flu ffy and frac- 
tal ag g regates are conipacted upon collision (iDomin ik & Tielen^ 
1 19971: lBlum& WurmlboOOh . A further increase in the impact 
velocity leads to compaction. However, eventually the kinetic 
energy is high enough to break contacts between individual 
monomers. Erosion then starts to remove parts of colliding ag- 
gregates. As the relative velocity increases further, erosion also 
becomes stronger, ultimately leading to the destruction of the 
aggregates. 
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Fragmentation of aggregates is a major obstacle in planet for- 
mation theory. Dust cannot grow all the way to form planetesi- 
mals because it gets destroyed once relative velocities become 
violent enough to disrupt aggregates. Similarly, radial drift may 
prevent growth by removing particles from the disk once they 
gro w to a certain size and spiral towards the central star. 

iJohansen et al.l (|2006') show that, in the presence of turbu- 
lence, planetesimals might be produced by gravitational collapse 
of clumps of meter sized boulders. These clumps are generated 
by high-pressure turbulent eddies that trap and concentrate parti- 
cles, which then become gravitationally bound. Before this pro- 
cess can take place, particles must already have grown by 18 
orders of magnitude in mass. The only feasible way to do that is 
by collisional sticking. Because collisional fragmentation of ag- 
gregates may prevent growth of large, meter sized aggregates, it 
is crucial to fully understand mechanisms involved in collisions 
of porous agg regates. 

Recently, iBrauer et"an (l2008l) have studied the growth of 
dust particles in a local density maximum caused by the evap- 
oration front at the snow line. This local pressure maximum in 
the disk can accumulate dust particles and reduce relative veloci- 
ties considera bly. In the midplane of the disk, weak er turbulence 
Cieslall2007h results in lower relative velocities. iBrauer et al.l 
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2008h have shown that for the fragmentation threshold veloc- 



ity of at least 5 m/s, the growth can proceed to large boulders of 
up to several 100 meters in size. However, a more realistic value 
of the threshold vel ocity, for aggregates made of micron- sized 
grains, is ~ 1 m/s (' Blum & Muenchlll993l: iPoppe et al.l I2OOOI: 
iLangkowski et al. 2008^^ 

[Blum & WurmI (QOO(f} performed laboratory experiments of 
collisions of dust aggregates. They studied impacts at a wide 
range of energies, from low energy (perfect sticking), through 
restructuring to fragmentation of microscopic aggregates. Their 
results with respect to res tructuring threshold are in agreement 
with theoretical findings of lDominik & Tielensl(ll997l) . The frag- 
mentation energy, however, diff'ers, which is a conseque nce of 
the discrepancy in sticking velocity. iPoppe et aP (l2000l) mea- 
sured the sticking veloci ty of micron-sized silica grains and 
found 1.2 m/s. Chokshi et al.l (1 19931) and Dominik & Tieleni 
(1 19971) on the other hand derived theoretically a much lower 
velocity, in consistent with experimen ts. This problem was ad- 
dressed by iPaszun & DominilJ (l2008l) . who included an addi- 
tional energy dissipation channel in order to match the experi- 
me ntal results. For details se e section [2l 

iBlum & MuenchI (1 19931) studied collisions of macroscopic, 
mm-sized, aggregates at velocities between 1 and a few m/s. In 
this case, however, a diff'erent behavior was observed. Although 
particles were very porous (up to 4% filling factor) they did 
not observe restructuring. Instead, aggregates were bouncing off' 
each other or, for faster impacts, fragmenting. 

Very energetic collisions between large mm-s ized a nd cm- 
sized aggregates were studied by Wurm et al. (l2005l) . They 
showed that the fragmentation observed at high velocity impacts 
turns into a net growth of 50% at velocities above 13 m/s. The 
distribution of fragments at velocities of about 20 m/s followed 
a power-law with a slope of -5.6 ± 0.2 for larger fragments and 
wa s ff at for the smallest e j ecta. 

iFujiwara et al.l (1 19771) experimented with solid basalt rocks. 
High velocity impacts of a few km/s result in a power-law dis- 
tribution of small fragments n(m) ~ m^. The slope of the distri- 
bution was found to be ^ = -1.83. They distinguished several 
collisional outputs, depending on the target size: 



ISironol (|2004|) has developed a smoothed particles hydrody- 
namics (SPH) model to simulate meter sized and larger aggre- 
gates. A single particle in this model corresponds to a porous 
material described by compressive strength, tensile strength, 
den sity, and sound speed. This method was also used by 
iSchafer et al.l (l2007h . 

Although collisions of dust particles have been studied ex- 
perimentally and theoretically, no one has formulated a quan- 
titative recipe describing both mass distribution and structural 
properties of the collisional output, using a model based on em- 
pirical results. Both the distribution of masses and the compact- 
ness of fragments is required to fully understand the growth of 
dust. Here we present an extensive parameter study of many col- 
lisions of small dust aggregates. We provide a recipe for mass 
distribution and compactness of fragments. 

In Sect. [21 we brieffy present the model we adopt to simulate 
collisions of dust aggregates, it's strengths and limitations. In 
Sect. [3] we present findings of our study and discuss them quali- 
tatively. Section IHprovides quantitative description of our results 
in a form of a collision recipe. We end this paper with interesting 
conclusions in Sect. \5\ 

2. The Model 

The simulations presented in this work are done using the N- 
body dynamics code SANEQ. Our model treats all monomers 
(also referred to as grains or particles) in the agglomerates (also 
referred to as aggregates, clusters, or particles) individually. 
Since we are currently not interested in long range interactions, 
electrostatic, magnetic and gravitational for ces are not included 
even though the code can handle them (iDominik & Niiboldl 
[2002h . We calculate motion of individual monomers that inter- 
act with each other via attractive van der Waals surface forces 
(iJohnson et al.ll 19711) . 

The presence of the attractive surface on deformable par- 
ticles inevitably leads to several energy dissipation mecha- 
nisms. The particles, when in contact, may roll over each other. 
This rolling motion is co untered by a rolling friction force 
(IDominik & TielensI [T995b . causing an energy loss. The same 
happens in the case of a sliding motion. The contacts may 
shift which again is work d one against the sliding friction force 
(IDominik & TielensI Il996h . Beside that, every time a contact 
between two mono mers is broken, the elastic energy stored 
in it is partially lost (iJohnson et al.l 1 19711: IChokshi et all 1 19931: 
IDominik & TielensI 1 19971) . Some energy m ay also be lost due 
to a t wisting motion of particles in contact (IDominik & TielensI 
[19971) . 

We also include an additional energy dissipation channe l 
in order to fit the experimental results by iPoppe et aP (l2000l) . 
This process increases the sticking velocity from about 10 cm s"^ 
(Chokshi et al. 1993; Dominik & Tielenslll997h to about 1 m s'^ 
(iPoppe et al.ii200Q) . Since the meas ured attractive f orce a grees 
well with the t heory provided by IJohnson et al.l (Il97ll) and 
iDerjaguin et al.l (Il975l) . the diff'erence in sticking velocity points 
to additional energy losses in collisions. While our model does 
not specify what this mechanism is, a candidate would be plas- 
tic deformation of surface asperities on nm scales. In order to 
achieve agreement with the experimental results, a mechanism 
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dissipating the energy upon the first co ntact of two particles 
was introduced (iPaszun & Dominil3l2008l) . For details regarding 
the im ple mentation, we ref er the rea der to Ipominik & Niiboldj 
(l2002h and lPaszun & Dom inik (20 08|). 

Although iPaszun &"Dominik (120081) tested the model ex- 
tensively against the laboratory experiments and found a good 
agreement, this model has limitations as discussed below. 

2.1. Limitations 

Although our model presents a new approach to study collision 
dynamics of aggregates, the following limitations apply. 

- Aggregate size 

Two monomers in contact oscillate in relative distance due 
to the competition between attractive forces that hold them 
together, and the elastic force that pushes them apart. For 
micron- sized grains (the monomers we consider here have 
radii of 0.5yum) the vibration frequency is on the order of 
GHz. A correct simulation of this dynamical system requires 
that we must resolve the shortest timescales. These very 
short simulation time steps limit the number of monomers 
in the system we can model. The largest aggregates we can 
simulate are made of 10^ grains and take several hundreds 
of hours of a CPU time. In this study we model aggregates 
made of up to 1000 monomers. 

The collisional outcome depends foremost on the strength 
of individual contacts in the impacting aggregates, thus the 
properties of monomers (size, composition). The aggregate 
size determines the amount of energy that can be dissipated 
(£^br times the number of monomers). Note, however, that 
laboratory experiments of much larger aggregates (about 100 
micron in size), made of billions of monomers, show phe- 
nomena n ot seen in the smallest a ggregates. An example is 
bouncing (iLangkowski et aD [2008h . Although this is an im- 
portant eff'ect that significantly aff'ects the growth of dust ag- 
gregates, it is still poorly understood. 

- Irregular grains 

Our model assumes spherical monomers, as the simulation 
of irregular grains is computationally very expensive and, 
therefore, is impractical. Moreover, we can directly com- 
pare our findings to the results of laboratory experiments 
um & Wu rm 2000). Here we briefly introduce the possi- 
ble eff'ect irregular grains may have on the collision outcome. 
The strength of an aggregate strongly depends on the size 
of the contact between monomers. In the case of randomly 
shaped grains, contact is established between surface asper- 
ities and is very much limited. This reduces the strength of 
an aggregate and leads to disruption at lower impact ener- 
gies. Moreover, irregular monomers can form more than one 
contact with each other, which means that individual con- 
tact points must break before restructuring can occur. In this 
way, chains of irregular particles can be more rigid than those 
made of sph e rical p articles. 

iPoppe et aP (l2000l) performed laboratory experiments on 
both spherical and irregular grains. They showed that irregu- 
lar particles can stick to a flat surface at much higher veloci- 
ties than spherical grains (1.5 to 2.3 m s"^ for spherical silica 
grains and 5 to 25 m s"^ for irregular enstatite monomers). 
This suggests that additional energy dissipation may occur 
during a collision. The mechanism, however, is still unclear. 
As the physical processes involved in collisions of irregular 
grains are unknown, we leave them for a future investiga- 
tion. However, experiments have shown that the overall ef- 



fects ofl' collisions are well reprodu ced between experiment s 
with spherical or round monomers (ILangkowski et al.ll2008l) . 
Therefore, we believe that the general conclusions of our 
study here will also hold for less perfect monomers. 

2.2. Setup 

A few examples of our aggregates, made of 1000 monomers 
each, are presented in Fig.[T] These particles consist of equal 
size monomers. These aggregates are constructed usin g two dif- 
ferent techniques. The first method, developed by Fili ppov et al.l 
(|2000|) . allows to create fractal aggregates of any specified frac- 
tal dimension Df. This sequential tunable particle -cluster ag- 
gregation method forms agglomerates by successive addition of 
identical spherical particles. The fractal aggregates made ac- 
cording to this algo rithm strictly obey the fractal scaling law 
(IFilippov et al.ll2000l) 

where is the radius of gyration and kf is the fractal pre-factor 
that is related to the central packing density. For the aggregates 
used in our study, we choose a pre-factor = 1.6 and fractal di- 
mensions of Df = 1.5, Df = 2.0, and Df = 2.5, respectively. We 
measure the filling factor of aggregates by considering a sphere 
with radius Rg-, having a cross section nR^. which is equal to the 
projected surface of the aggregate, averaged over all angles (see 
section l234l) . We then compute the filling factor (po- as the ra- 
tio of the solid volume in the aggregate and the volume of the 
equivalent sphere. With this definition, aggregates of the men- 
tioned fractal dimensions, made of 1000 monomers, have filling 
factors (po- = 0.07, (po- = 0.09, and (pa- = 0.16, respectively. 

The second method we use to produce aggregates is particle 
cluster aggregation (PC A). We successively add monomers from 
random directions. This method produces aggregates that, in the 
limit of very large sizes, have a filling factor of 0^- = 0.15. An 
aggregate made of 1000 monomers has, however, a lower filling 
factor of (po- = 0.13. This is an eff'ect of high porosity of the sur- 
face layer. In very large aggregates, this re gion will be negligib ly 
thin compared to the size of an aggregate (iPaszun et al .1120091) . 



2.3. Parameter space and expected scaling 

To provide a qualitative and quantitative description of aggre- 
gate collisions, we explore an extensive parameter space. This 
provides an insight into the eff'ects of diff'erent parameters on the 
outcome of a collision. As the final purpose is to provide a recipe 
for a collision between two aggregates, we limit the range of our 
parameters to realistic values. 

The relevant parameters influencing the outcome of a colli- 
sion are: 

- collision energy, 

- pre-collision compactness of the aggregate, 

- mass ratio of colliding particles, 

- impact parameter, 

- material properties of monomers. 

The values used for the difl'erent parameters are shown in Tab. [T] 
The columns in the table refer to (l)-relative collision velocity 
V, (2)-pre-collision filling factor of an aggregate (po-, (3)-impact 
parameter 5 b/bjnax, and (4)-mass ratio of colliding aggregates 
mi/m2. Although not all combinations are simulated, each pa- 
rameter set used it computed for six diff'erent orientations in or- 
der to obtain a fair average of possible collision outputs. 
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Fig. 1: Sample aggregates used in our parameter study. Each of particles is made of 1000 monomers. Annotations show ( 
circles define the outer radius /?out^ while gray circles indicate the projected surface equivalent radius R^^ (see text). 



^for each particle. Black 



Table 1: 

study. 



Impact parameters and aggregate properties explored in this 



V [m/s] 


(per 




mi/m2 


(1) 


(2) 


(3) 


(4) 


0.01 


0.07 


0.0 


10-^ 


0.05 


0.09 


0.25 


5 • 10-2 


0.1 


0.122 


0.5 


1 • 10-1 


0.2 


0.13 


0.75 


2-10-1 


0.3 


0.155 


0.875 


4-10-1 


0.5 


0.16 


0.95 


6-10-1 


0.75 


0.189 




8 - 10-1 


1.0 


0.251 




1 


2.0 








4.0 








6.0 








8.0 








10.0 









2.3.1 . Collision energy 

The basic effect of impact energy is intuitive: The more en- 
er gy is provided to th e syste m, the m ore violent the out come 
is. iDominik & TielensI (Il997h and later IWada etaP (l2007h pro- 
vide a simple recipe of collisional output as a function of en- 
ergy. Both these studies were limited to two dimensional ag- 
gregates. Their recipe predicts energy thresholds for processes 
such as erosion, compression, and fragmentation. To understand 
them, one needs to define the rolling energy and the breaking 
energy of a contact. The first one is the energy needed to roll 
two monomers over each other by 90 degrees, and it represents 
en ergy the related to restruct uring of an aggregate. It is defined 
as (iDominik & TielensI! 1997h 



^roll 



6n yR^cnu 



(2) 



where y is a surface energy, R is reduced radius of the two 
monomers in consideration, and ^crit is a critical displacement at 
which the rolling becomes irreversible and energy is dissipated. 
From measurements of the rolling friction between spheres it 



was shown that this quantity should be se t to about 20A in o rder 
to reproduce the measured energy losses (iHeim et al.lll999h . 

The second important quantity is the energy needed to sepa- 
rate two connected monomers and is defined as 



^br — 1.8/^c^c 



(3) 



Here Fc is the pull-oflT force, minimum force needed to discon- 
nect two monomers and it is 



F, = 37iyR 



(4) 



where is the critical displacement from the equilibrium posi- 
tion (Dominik & Tielens 1997). Thus, two monomers in contact, 
when pulled apart, will separate at this position, is defined as 



1 < 



2 61/3/?' 

with the equilibrium contact radius 

/97rr^^i/3 



ao 



(5) 



(6) 



Putting these equations together we see that the critical energy 
is given by 

5/3^4/3 

with the dimensio nless constant A = 4 3. 

The recipe of IDominik & TielensI (1 19971) is confirmed ex- 
perim entally with one important modification (iBlum & WurmI 
|2QQQ|) . The energy scaling should be applied according to values 
determined empirically (i.e. ^crit ~ 20A and £'br corresponding 
to the experimental results). In order to satisfy this requirement, 
we introduce a scaled version of £^br. In this case the constant 
A is higher such that the stic king threshold in our model is in 
agreement with experiments (iPoppe et al.ll2000l) . This energy is 
given by 



(8) 



The recipe by IDominik & TielensI (1 19971) is summarized in 
Tab. [21 Low energies are insuflScient to cause any visible re- 
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Table 2: The collision recipe from lDominik & TielensI (119971) for a 2D 
case. 



Energy 


Outcome of Collision 


-^impact ^ 5£'roll 
^impact ~ 5 £"1-011 

^impact ~ '^c-^roll 
-^impact ~ 3/tc-£'br 

^impact ^ lO'^c^br 


Sticking without restructuring 
Onset of restructuring local 
to the impact area 
Maximum compression 
Onset of erosion (start to 
lose monomers) 
Catastrophic disruption 



structuring. Before any restructuring occurs, contacts between 
monomers are displaced elastically. Only when the critical dis- 
placement is exceeded, an irreversible motion (i.e. rolling or 
sliding) follows. When the energy is sufficient to roll several 
contacts by a significant angle (-90 degrees), monomers in the 
impact region begin to roll and visible restructuring occurs. 
The amount of restructuring increases with increasing energy. 
Maximum compression is typically reached, when all monomers 
have enough energy to roll by 90 degrees. 

Monomers begin to be removed of an aggregate when the 
energy per contact reaches 3 times the breaking energy. A large 
fraction of the energy is then dissipated by rolling, but the excess 
is used to remove a few grains. As the impact energy increases, 
the amount of erosion increases and becomes catastrophic when 
each contact can access over 10 times the breaking energy. 

Although this recipe is reasonable, it may diff'er when par- 
ticles collide at an impact parameter other than Z? = or when 
the mass ratio of the two col lision partners is no t 1. Th e latter 
case was also considered by iDominik & TielensI (1 19971) . They 
provided a recipe for collisions of a small grain with a cluster. 

In this study we also explore the energy range from the hit- 
and- stick regime up to catastrophic destruction. This range of 
energies can be found by looking at a variety of astrophysical en- 
vironments. For example, in protoplanetary disks, small particles 
are well coupled to the gas and have rather low relative veloci- 
ties. In the cores of molecular clouds, however, low gas densi- 
ties make it possible for particles to decouple from the gas even 
at small sizes. In this case, turbulent gas motions may lead to 
collis ions velocities beyond the fragmentation limit (lOrmel et al.l 
I2009h . 

2.3.2. Mass ratio 

The mass ratio is sampled within the available range of our 
model (10"^ < mi/m2 < 1). As the largest aggregate simulated 
in this study is made of 1000 monomers, the lowest mass ratio 
we consider is mi /m2 = 10"^ (a collision of a big aggregate with 
a monomer). Although we simulate collisions between particles 
with diff'erent mass ratios, we present only two limiting cases to 
illustrate the importance of the mass ratio. 

This parameter influences the energy distribution during a 
collision. In the case of equal-mass impactors, the energy is 
spread over most of the monomers, causing a global efl'ect. 
However, small projectiles act locally and concentrate the en- 
ergy into a small volume close to the impact site. Thus, restruc- 
turing will be local and erosion can be expected at relatively low 
energies. 



2.3.3. Impact parameter 

Our sampling of the impact parameter b covers the range from a 
central impact up to a grazing collision, where the impact param- 
eter b equals 95% of the sum of the outer radii of two colliding 
aggregates 7?out,i + ^out,2. We do not consider impact parameters 
above 95% as the irregular shape of aggregates generally causes 
such "collisions" to miss. 

The impact parameter is important as it can significantly 
change the outcome of a collision. In the case of a central im- 
pact, aggregates are pushed towards each other and are com- 
pressed. Collisions with large impact parameters on the other 
hand cause aggregates to connect in the outer regions only and 
result in stretching of aggregates as they move apart again. In 
this way, large impact parameters tend to cause tensile forces 
acting on the aggregates, while central impacts are dominated 
by compressive stresses. 

The importance of the off'-center collisions should not be 
underestimated. Due to the geometrical arrangement, impacts 
with large impact parameter are considerably more frequent than 
head-on collisions. In impact parameter averaging, grazing im- 
pacts have a lot of weight. 

2.3.4. Compactness parameter 

To describe the structural changes of aggregates we define a 
compactness parameter 0^- as 

where N is the number of monomers in the aggregate, tq is a 
monomer radius, and Ro- is the projected surface equivalent ra- 
dius (cf. Fig.[T]), defined as 

Rcr = j^. (10) 

with (T being the projected surface averaged over many orienta- 
tions. The inverse o f this filling factor was introduced earlier by 
lOrmel et aP (l2007l) as the enlargement factor i//. Note that this 
parameter alone is insufficient to fully describe the structure of 
an aggregate. An additional quantity is required to avoid ambi- 
guity, e.g. the fractal dimension. However, for small aggregates 
such as the ones in this study, the possibilities to distribute mass 
within the particle are limited. 

We sample the compactness parameter well within the ap- 
plicable range. Random close packing (RCP) of spheres pro- 
duces aggregates with t he filling factor of abouQ (po- ~ 0.635 
(IOnoda&Linigejll990l) that is the densest form one can ex- 
pec t, given the assumption of spherical monome rs. Ho wever, 
Blu m etal.1 ^006) and later Paszun & Dominik (l2008h have 
shown that aggregates being compressed can reach a maximum 
filling factor of about cpo- = 0.33. Higher compaction cannot be 
achieved in a static experiment of uni-axial compression. As up- 
per limit for the filling factor we use a slightly lower value of 
(pa- = 0.25. As lower limit, on the other hand, we use aggre- 
gates formed in the Brownian growth phase, where in the pres- 
ence of rotation aggr egates with the fractal dimension of about 
Df = 1.5 are formed ause &Blumll2QQ4l: IP aszun & DominikI 
'2006'). Our largest aggregate of fractal dimension 1 .5 has a filling 
factor 0^ = 0.07. 

^ Although the value of 0.635 correspond to filling factor defined in 
respect to the outer radius Rout and not Rcr, both radii are equal for these 
compact aggregates. 
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0.01 I I '■■ , , , ^ , , , I 

1 10 100 

Rout/Ra 

Fig. 2: The filling factor as a function of the ratio Rout/Ra of the outer 
radius of an aggregate over its projected surface equivalent radius. Each 
point corresponds to an aggregate (see text) of different structure and 
mass and thus different ratio of the two radii. The dashed-dotted line 
shows the least square fit of the power-law to this data. A sharp transi- 
tion between compact and fluffy particles occurs at about i?out/^cr = 1.2 
(dashed line). Compact aggregates have a smaller ratio of the two radii. 
Their filling factor is approximated by a separate power-law (dotted 
line). 



It is important to note the discrepancy between the radii used 
in the definition of the impact parameter space and the filling 
factor (po-- The impact parameter is defined in terms of the outer 
radius 7?out, that is the radius of a sphere enclosing an entire ag- 
gregate and centered in its center of mass. The filling factor, 
however, uses the projected surface-equivalent radius Ro-. For 
compact aggregates, these two radii are very similar, while the 
outer radius Rout becomes higher as the filling factor of an ag- 
gregate decreases. For completeness, we empirically determine 
the relation between the two radii. Figure [2] shows how the ra- 
tio Rout/R(T influences the filling factor (pa-- To obtain this rela- 
tion we constructed (using Eq. ([T])) many aggregates of differ- 
ent structures (Z)f, kf) and masses (N). Interestingly all data is 
very well confined along a simple curve. As the filling factor 
seems to depend on the mass of aggregate as 0o- ^ 
plot (pa-N^'^^ in Fig.[2]to collapse all data on a single curve. This 
mass dependence is further discussed in Sect. 13. 2[ Compact par- 
ticles in Fig. [2l(7?out/^£r < 1 .2) have a similar outer radii Rout and 
projected surface equivalent radii Rg-. They show little depen- 
dence of the Rout/R(r ratio on the filling factor. Fluffy aggregates 
(Rout/R(r > 1.2) on the other hand show a power-law relation of 
the filling factor on the ratio of the two radii with a slope of -0.3. 
The fitted power-law is presented in Fig.[2las dashed-dotted line. 
The complete relation that holds for ratio Rout/Ro- > 1.2 is given 
as 

0<,= i.2i(^pyv-o-3l (11) 

2.3.5. Properties of monomers 

The composition of monomers and their size strongly affect the 
strength of an aggregate. Both these parameters determine the 
breaking energy of two grains in contact and thus regulate the 
energy dissipation during a collision. Since normalization of the 
impact energy by the breaking energy puts monomer properties 
out of the equation, we study here only one monomer size and 



material. The physical parameters of Quartz used in our model 
are presented in Tab. [3] Other materials will be the subject of a 
future study. Our monomers are silica spheres with a diameter of 



Table 3: Properties of monomers used in this study. Material properties 
correspond to Quartz monomers or radius tq = 6 • 10~^ cm. 

y [erg/cm^] [dyn/cm^] ^crit [cm] p [g/cm^] 
25 2.78-10^^ 2 • 10"^ 165 



1.2 yum. The elasticity modulus E* is defined as 

where Ei and v/ are the Young's modulus and the Poisson ratio 
of i-th monomer, respectively. 

3. Results and discussion 

In this section we present results of our parameter study. We 
describe the collisional outcome in terms of the fragment mass 
distribution - growth versus fragmentation. Moreover, we keep 
track of the structure of the fragments formed in such a colli- 
sion and present the collisional evolution of the structure of dust 
aggregates. 

Our study spans a wide range of parameters. Here we select 
two specific cases for a detailed discussion of the effects seen 
in our stud}0. Collisions of compact aggregates (filling factor of 
(pa- = 0.251 for smaller aggregates and (po- = 0.16 for larger par- 
ticles) are presented and compared with fluffy, fractal aggregates 
(filling factor of (po- = 0.155 for small aggregates and (po- = 0.09 
for bigger ones). We use these two cases to illustrate the influ- 
ence of compactness on the collisional outcome. 

3.1. Fragment distribution 

The products of a collision are generally quantified in terms of 
the fragment mass distribution. In the case of sticking, the result- 
ing mass distribution contains one single element, with the mass 
given by the sum of impactor and projectile. With the onset of 
erosion, a second component appears - the distribution of small 
fragments, usually represented as a power-law of particle mass. 
While these are initially two clearly separated components, they 
can connect in increasingly destructive collisions. 

3.1 .1 . The effect of impact energy 

Basic mass distribution components Although the impact en- 
ergy seems to be the main quantity setting the mass spectrum af- 
ter the collision, the collisional outcome depends very much on 
all parameters presented in Sect. 12. 3[ The largest collisional rem- 
nant is presented as a function of the impact energy in Fig. [3] At 
low collision energies, aggregates stick perfectly and the largest 
fragment contains the mass of both colliding particles. An in- 
creasing impact energy causes onset of erosion at the energy of 
about E = 0.1NE\)Y. The fragmentation occurs at the energy of 
about E ~ NE\)Y, and depends also on the internal structure of 
the colliding aggregates and the impact parameter. 

^ The quantitative recipe presented in Sect. lU does of course make 
use of the entire parameter study. 
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Fig. 3: The mass of the largest collision remnant for central and offset 
collisions as a function of the impact energy. Diamonds correspond to 
compact aggregates, while triangles correspond to fluffy particles. Error 
bars indicate standard deviation and each point represents a mean value 
from 6 individual simulations (see Sect. l2.2l ). 



Compact aggregates with densely packed monomers = 
0.251) can sustain higher energies regardless of the impact pa- 
rameter. Grains packed close to each other undergo more interac- 
tion, resulting in internal energy dissipation that is more efficient 
than in the case of loosely packed monomers ((po- = 0.155). 

Interestingly, erosion occurs in off-center collisions already 
much earlier than in central collisions, while shattering has the 
opposite behavior: It occurs more readily in central collisions. 
Central impacts can cause catastrophic destruction at energies 
above E = 10A^£^br, while off-center collisions at this same en- 
ergy usually cause two large fragments to remain and therefore 
must be classified as be classified as erosion. Both effects are 
caused by the distribution of the available impact energy within 
the colliding aggregates. In off-center collisions, the energy is 
concentrated into a small region of the aggregates and does break 
grain-grain connections locally near the point of impact. The re- 
maining kinetic energy is then carried away by two massive frag- 
ments. In the limiting case of (b = /?max) the interaction of the 
two aggregates occurs between only two monomers, if at all. 

The degree of fragmentation during off-center impacts de- 
pends then on the packing density of grains in the two aggre- 
gates. Porous aggregates (cpo- = 0.155) show that even at very 
high energies, the energy is not efficiently absorbed by the ag- 
gregates, saving them from being shattered. However, in the case 
of compact aggregates, the mass of the largest collision remnant 
is significantly higher than in the case of central impact, but it 
decreases with an increasing impact energy. 

The full mass distribution A more complete picture is pre- 
sented by studying the full distribution of fragment masses. 
This is illustrated for central impacts in Fig.lH Collisions at 
intermediate speeds (top panels) erode aggregates and produce a 
number of small fragments. Both the power-law and the largest 
fragment components are present. Depending on the degree of 
erosion they can be separated or connected. However, the slope 
of the power-law remains more or less unchanged. 
The distribution of small particles is defined as 




f(m/mo)(m/mo) oc (m/mo)^ 



(13) 



10^ 10^ 
m/rrio 



Fig. 4: Mass distribution resulting from head-on collisions of A&C - 
fluffy (0^ = 0.155) and B&D - compact (0^^ = 0.251) aggregates at 
four different impact energies. 



with the slope q depending on several parameters, including the 
internal structure of aggregates and the impact energy. Initially, 
too few small particles are produced and the power-law cannot 
be determined (e.g., for off-set impacts see Fig.|7^&b). In such a 
case we assume that the distribution is flat with the slope q = 0, 
meaning equal mass per logarithmic mass interval. As the ero- 
sion progresses with increasing impact energy, the fitted slope 
remains at an approximately constant low value. Small particles 
slightly dominate the mass within this power-law distribution, 
since the slope is about q ^ -0.3 in the case of fluffy aggregates 
(Fig.|4^), and about q ^ -1.2 for compact aggregates (Fig. St). 
These slopes remain constant up to energies of about E = NE]^^. 
Faster impacts provide sufficient energy to break all contacts in 
an aggregate. Therefore, full fragmentation sets in, which can 
significantly alter the fragment distribution. The head-on colli- 
sion at energy ofE = 3.66A/^£'br completely shatters fluffy ag- 
gregates with the largest fragment having only about 15% of the 
total mass. The large fragment component disappears altogether, 
shifting the entire mass into the power-law component with the 
now steeper slope q ^ -I. Further increase in the energy re- 
sults in heavier damage and a still steeper slope. At the energy 
ofE= 14.6A/^£'br the slope reaches q ^ -2.5 (see Fig. [4]:). The 
steepening slope can be interpreted as monomers and very small 
fragments in the fragment distribution becoming dominant. 

Although a similar trend is observed for compact particles 
((po- = 0.251), the fragmentation is not as effective as in the case 
of aggregates with 0^- = 0.155 (see Fig.|4]l). A head-on impact 
of compact particles at the energy of E = 3.66A/^£'br can still 
be classified as erosion. The large fragment component is still 
present and contains a significant fraction of the mass. This com- 
ponent is broader and is connected with the power-law distribu- 
tion of small fragments, an effect which affects the determination 
of the slope. The slope appears to decrease to about q ^ -0.3. 
However, an increase in the collision energy to E = 14.6A/^£^br 
shatters the aggregates, leaving only the power-law component 
with a steeper slope of about ^ -1.7 (see Fig.|4jl). 

The relation of the energy and the slope of the power-law dis- 
tribution in central collisions is presented in Fig. \5\ The left panel 
(a) shows the results for aggregates with compactness parame- 
ter (po- = 0.155. The initially shallow and constant slope begins 
to steepen once the impact energy increases above E = A/^£^br- 
This behavior is similar for aggregates of different masses or for 
different mass ratio impacts. Note that for very weak erosion. 
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Fig. 5: Slopes of the power-law size distribution of small fragments 
as a function of impact energy. Different symbols correspond to dif- 
ferent masses of the colliding aggregates. The left panel (A) corre- 
sponds to fluffy aggregates (0^^ = 0.09 . . . 0.155), while the right panel 
(B) shows the result of collisions between compact aggregates (0^^ = 
0.16... 0.251). 
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Fig. 6: Mass distribution resulting from central collisions between a 
monomer and an aggregate made of 1000 grains. The presented dis- 
tributions result from collisions with compact aggregate with the filling 
factor of = 0.16. Note that the energy is normalized to the reduced 
number of monomers = MA^2/(M + A^2)- 



the slope of the power-law cannot be determined and a value 
of ^ = is assumed. Shattering becomes catastrophic when the 
impact energy is increased by an order of magnitude. The slope 
steepens beyond q = -I, meaning that small particles dominate 
the mass spectrum. For small aggregates the distribution may be 
as steep as ^ = -2.5. 

Compact particles behave in a similar manner (see Fig.[5]3). 
In this case, however, the steepening occurs at slightly higher en- 
ergies of about E = 5NE\yY. Moreover, the steepening is limited 
in the explored energy range, and reaches values of ^ = -1.7. In 
this case, low energies also result in shallow slopes of the distri- 
bution and are again assumed to be ^ = when the data was too 
scarce to make a fit. 



3.1 .2. The effect of impactor-to-target mass ratio 

Below we discuss the effect of the mass ratio on the collision 
outcome. Although intermediate energies are required to erode 
particles in collisions of equal mass aggregates, this is not the 
case for impacts with high mass ratios, where the impact energy 
is localized to a small region, leading to erosion even in slow 
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Fig. 7: Mass distribution as resulted from off-center (b = 0.75 Z?max) 
coUisions of A&C - fluffy (0^ = 0.155) and B&D - compact = 
0.251) aggregates at four different impact energies. 



collisions. Figure [6] shows the fragment distribution produced in 
collision of a monomer with an aggregate made of 10^ grains. In 
these cases, the energy is sufficient to break only a low number 
of contactfl and results in erosion. That same energy applied in 
a collision of equal mass aggregates results in a perfect sticking 
without any mass loss (cf. Fig. [3]). This difference is a conse- 
quence of very localized energy input. The small particle (in this 
case a monomer) carries sufficient energy to break a number of 
contacts. This energy is transmitted locally to a limited number 
of grains rather than distributed over the entire target aggregate. 
The resulting ejecta can, due to the small physical size of the 
projectile, easily escape. 

Similarly, off-center collisions (see Fig. [7]) are characterized 
by distributions that resemble the erosion case. A highly pro- 
nounced large fragment component coexists with a power-law 
distribution of small fragments. The slope of the power-law is 
independent of the impact energy. It remains at about q ^ -1 
for fluffy aggregates and q ^ -0.6 for compact particles even at 
very high energies above E = 10A/^£'br. At low impact energies, 
the resulting distribution are dominated by the largest fragment 
component and the slope of the power-law component is unde- 
termined due to the low number of fragments. 

Fluffy aggregates, when colliding at an energy of £" = 
0.23A/^£'br and above, produce two large fragments that con- 
tribute to the big particles component and a few small parti- 
cles. Compact aggregates can still grow at an energy of £" = 
0.23A/^£^br, producing one large fragment containing most of the 
mass. At energies of about E ^ NEi^^ the erosion increases and 
the growth is stopped. The large fragment component contains 
now two large remnants of the initial aggregates. An increase 
of the collision energy causes mainly an increase in the mass 
in the power-law component, while the largest fragment remains 
weakly affected. In the case of compact aggregates ((po- = 0.25 1), 
only very high energy impacts can significantly reduce the mass 
of the largest fragments (see Fig. [3 and Fig.|7]l). For aggre- 
gates with lower packing density (0^- = 0.155), the mass of 
the largest collision remnant seems to stabilize at a value of 
about mmax ~ 0.5 (mi + m2) once the impact energy exceeds 



^ Note that in Fig. [6l the energy is normalized to the reduced number 
Therefore, the energy per contact is a factor 



of monomers = 



of 10^ lower than in the case of a collision of particles of equal mass. 
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Fig. 8: Schematic plot of the evolution of the fragment distribution as a 
function of different parameters. 




Fig. 11: The geometrical filling factor (per as a function of impact energy 
for head-on (solid lines) and off'-center (dashed lines) collisions. Black 
lines correspond to aggregates with lower initial filling factor of = 
0.155, while gray lines correspond to higher filling factor of 0o- = 0.251. 



E = NE\)Y. Very high energies influence the fragment distribu- 
tion very weakly (see Fig.[3]and Fig. [7]:). 

The main efl'ect of the impact parameter is that the energy 
is not transported very efficiently into aggregates. In the case 
of a central collision the kinetic energy is naturally transported 
efficiently into both aggregates, as the interaction spreads from 
the center (i.e., region right in between the aggregates) outwards. 
The monomers are pushed closer together, actively taking part 
in the energy dissipation. An increase in the impact parameter 
results in a decrease in the number of actively interacting grains. 
This means that fewer grains actually collide resulting in less 
fragmentation. 

Our results show that the main factor determining the colli- 
sional outcome is the impact energy E. Other parameters influ- 
ence the way this energy is transported to and distributed over 
the available monomers. Figure [8] shows a schematic picture of 
how the fragment distribution changes with variation of impact 
energy E, compactness parameter (pa-, impact parameter b, and 
the mass ratio mi /m2. The arrows indicate schematically how the 
position and scaling of the difl'erent components of the fragment 
mass distributions shift as parameters vary. 

3.2. Aggregate structure 
3.2.1 . Central collisions 

The outcome of a collision can also be quantified in terms of 
the internal structure of the resulting aggregates. This quantity 
is very important as it determines the aerodynamic properties 
of aggregates and thus their relative velocities. Figure [9] shows 
examples of the effect of collisions involving aggregates with 
(po- = 0.155. Central collisions evolve from the hit and stick 
regime through compaction and flattening to fragmentation. As 
the impact parameter increases, this path changes. Then the hit 
and stick regime is followed by the stretching regime, where ag- 
gregates are pulled into an elongated shape, exposing more sur- 
face in the process. Eventually, suflftcient energy is reached to 
break fragile connection between the colliding aggregates, at the 
location of the highest stress. At this moment, two aggregates 
of similar mass are created. Depending on the impact parameter 



and energy, erosion may also occur in the vicinity of the interac- 
tion area. 

This picture depends also on the initial compactness of the 
colliding aggregates. Particles with open structures (Fig. [9]) are 
subject to drastic restructuring, while more compact ones are 
characterized by higher strength against compacting and ten- 
sile forces. For comparison we present the results for aggre- 
gates with geometrical filling factor of (pa- = 0.251 in Fig.fTOl 
In this case the produced aggregate is not stretched as much as 
in the case of more fluflTy aggregate. Additionally, the erosion is 
much stronger at higher impact parameters. This is caused by the 
closer packing of monomers that enables more eflftcient energy 
transfer into the aggregates. 

The quantitative picture is presented in Fig.[TT] The geomet- 
rical filling factor is plotted as a function of the impact energy for 
central and oflT-center collisions. All post-impact aggregates are 
in the hit and stick energy regime, below a few times the rolling 
energy Ero\\. This is consis tent with the description provided by 
iDominik & Tielensl (1 19971) . where the visible restructuring oc- 
curs at impact energies above 5 E^ow. However, the structure of 
the produced aggregate is dififerent than that of the initial parti- 
cles. Aggregates stick, forming dimer-like structure. Thus, they 
expose more surface area resulting in a decrease in the geometri- 
cal filling factor. This decrease depends on the initial structure of 
aggregates and appears to be more severe for compacter aggre- 
gates. In this case, particles cannot penetrate each other contrary 
to more open structure aggregates, where partial overlap hides 
some surface area. In fact, the diflference is rather small, on the 
order of 20%, for both filling factors presented in Fig.fTTl 

Further increase in the impact energy, beyond 5£'roii, re- 
sults in an increase of the filling factor. In the case of cen- 
tral collision aggregates undergo compression up to energy 
of about 0.5A^£'roii. At this energy, the maximum compaction 
is reached and this energy is close to threshold predicted by 
'Dominik & Tielens (1991) and confirmed experimentally by 
IBlum & Wurm (20001 The degree of compression depends on 
the initial compactness, as compact aggregates are more difl&cult 
to compress further. An aggregate with (per = 0.155 reaches its 
maximum compression of about (po- = OA 8+0.02, more compact 
than the pre-impact particles. A more compact aggregate with 
pre-impact filling factor (po- = 0.251 does not exceed it's ini- 




Fig. 9: Images of aggregates formed in collisions at different impact energies E and at different impact parameters b. The aggregates have an initial 
filling factor of = 0.155. Top panels show the initial setup. Lower panels present result of a collision. Images are in different scales to show 
different processes, i.e. compaction or fragmentation. 



INITIAL 



0.01 m/s 



1.0 m/s 



2.0 m/s 



8.0 m/s 




Fig. 10: Images of aggregates formed in collisions at different impact energies E and at different impact parameters b. Aggregates have the filling 
factor of = 0.251. Top panels show the initial setup. Lower panels present result of a collision. Images are in different scales to show different 
processes, i.e. compaction or fragmentation. 



12 



D. Paszun & C. Dominik: Collisional evolution of dust aggregates 



tial compactness after impact. The maximum compression ob- 
tained in our study is somewhat lower than obtained in a quasi- 
static compression (Blum et al. 2006; Paszun & Dominik 2008). 
This, however, can be an effect of small aggregate sizes, where 
the filling factor is strongly aff'ected by the porous outer layers 
(iPaszun et al.ll2009h . 

A further increase in the collision energy causes flattening 
and, formally, a decompaction of aggregates. Even thought the 
filled regions of the aggregates continue to show low poros- 
ity, the non-spherical global shape exposes more surface than a 
spherical structure would, leading to a decrease in the geometri- 
cal filling factor. In the case of flufl'y aggregates with 0^- = - 155, 
the flattening results in a small decrease in the compactness. At 
a collision energy of about N E^qw the maximum decompaction 
is reached and any further increase in the impact energy would 
lead to fragmentation. For more compact aggregates with initial 
0O- = 0.251, however, the flattening is stronger, as the filling fac- 
tor drops to about 0^- ~ 0-18. This filling factor can is similar to 
what we have seen earlier in flufl'y aggregates. Fragmentation is 
now only be reached at about 3 

The influence of the impact parameter 

Off-center collisions, on the other hand, show different results. 
Energies that lead to compaction in central impacts also cause 
restructuring at large impact parameters. In this case, however, 
particles are pulled apart and stretched. Thus, more surface area 
is exposed resulting in a strong decrease in the geometrical fill- 
ing factor (per- The stretching energy regime extends, however, to 
lower energies than the compression regime. Aggregates connect 
with a lower number of contacts that are pulled off. Therefore, 
lower energy is sufficient to disconnect the two colliding aggre- 
gates. The critical energy is about 0.25 E^qw, and may be slightly 
higher for compact aggregates. Above that energy, two parti- 
cles of similar mass are produced, accompanied by erosion. All 
small fragments produced both in central and off-center colli- 
sions closely follow a single power-law relation of 

0^ = {mlm^fl\ (14) 

The reason for this consistent behavior lies in the fact that the 
fragments are produced in highly dynamical events with energies 
close to breakup energies, allowing internal restructuring to lead 
to a scale-free structure. 

3.2.3. Schematic representation 

Figure [121 sketches a general picture of the structural evolution 
of aggregates. On the two axes, in arbitrary units, it shows filling 
factor and mass. We assume that initially, both impactor and tar- 
get have equal properties, located at the center of the plot where 
the solid black curve starts. Along that curve, the impact energy 
increases in steps, and the position of points on that curve shows 
possible structures of the largest post-impact aggregate. The first 
segment shows the effect of a hit-and-stick collision, to double 
mass and lower filling factor, independently of impact param- 
eter. At that point, the curve splits into two, for central (solid 
black) and grazing (gray) collisions. Higher impact energies can 
then, with equal mass, either increase or decrease the filling 
factor. Even higher energies move to the third point on each 
path. Grazing collisions return to the initial pre-impact proper- 
ties, while central collisions go through a maximum compres- 
sion point onto the fragmentation powerlaw (per = (m/mo)^.33. 

Small arrows in the sketch indicate how the different points 
shift around as model parameters are modified. 




mass 



Fig. 12: Sketch of the compactness evolution dependence on different 
parameters. See text. 

4. The Recipe 

In this section we put together the information regarding the col- 
lisional output in a form of a quantitative recipe. We provide 
the recipe in a form of tables that contain parameters needed 
to reconstruct appropriate distributions. Since the recipe is pro- 
vided as discrete parameters, the intermediate cases should be 
interpolated. Our recipe describes two limiting cases separately. 
Equal mass collisions affect aggregates globally. Therefore we 
refer to this sub-recipe as the global recipe. On the other hand 
we have high mass ratio impacts that result in localized changes. 
This sub-recipe is referred to as the local recipe. In the case of 
large mass ratio collisions, we used impacts of monomers onto 
aggregates composed of 1000 grains. In this case, only central 
collisions were used to predict also the outcome of offset im- 
pacts, because the effects are extremely similar. This assumes 
that the geometry of the impact of small projectile onto a larger 
target depends on a local surface. Thus central and offset colli- 
sions should be indistinguishable. The only necessary correction 
is to exclude missed collisions that are more likely at larger im- 
pact parameters. The distinction between the two sub-recipes is 
described later in Sect. 14.51 

4.1. Average over the impact parameter 

We present our recipe for aggregate collisions in an impact- 
parameter- averaged way. In this way, the recipe can be easily 
used in models that deal with size and structure distributions 
rather than individual particles. It is also possible to rewrite 
the recipe, keeping the impact parameter as one of the colli- 
sional parameters. This can be useful in a Monte-Carlo approach 
(lOrmel et al .1120091) where individual collisions are treated. 

Since we know the products of collisions at different offsets, 
we apply proper weights to the results. The weight for each im- 
pact parameter b is related to the fractional surface area of the 
ring with the radius b and width Ab. Thus, grazing collisions 
get the highest weight, while the central impact has the lowest 
weight. The quantity Q averaged over the impact parameter is 
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Table 4: Quantities provided by the recipe to reproduce the mass distri- 
bution. 

symbol description 

q The slope of the power-law component. 
Mr ratio of the mass in the power-law component 

to the mass in the Gaussian component. 
(Tq Width of the Gaussian component. 
Mq Mean mass of the Gaussian component. 



then given as 



<Q>b^ 



J^^"^'^ Q{h)27ih db 



(15) 



4.2. Hit and stick recipe 



Although our simulations include the hit-and-stick energy 
regime, this growth mechanism strongly depends on the mass ra- 
tio and should be treated separately. The hit-and-stick growth is 
of little importance when only an insignificant mass is added be- 
fore any restructuring oc curs. O therwise, an analytical prescrip- 
tion may be applied[ as inlOssenkopi ( 1993) : .Ormel et slL (2001. 
l2009h : lPaszun et al.l (l20Q9l) ^ 



4.3. Distribution of fragment masses 

The simulation results have shown that collisions produce two 
components. The first component - the power-law of small frag- 
ments - is produced in high energy, head-on impacts resulting 
in erosion and fragmentation. The second component - the dis- 
tribution of the largest fragments - is formed during low energy 
impacts, and also by grazing collisions at all energies. We model 
this component using (somewhat arbitrarily) a Gaussian distribu- 
tion of particle masses. The overall eff'ects seen in the previous 
chapter indicate that, as the energy increases, the Gaussian com- 
ponent will move its peak position to smaller masses, the power- 
law component will gain relative importance and will eventually 
steepen. 

The total set of parameters needed to reproduce such a dis- 
tribution is presented in Tab. IH The power-law component is de- 
termined by fitting a power-law to the first part of the distribu- 
tion, which contains small fragments. That power-law is then 
subtracted from the distribution. The power-law slope and mass 
ratio of the two components is known at this point. The remain- 
ing part of the distribution is then used to determine the last two 
quantities, namely the mean and the width of the Gaussian com- 
ponent. The mean mass Mq is simply the mass-weighted mean 
mass of fragments in the Gaussian component 



Mg = 



Jk 



log(mi+m2) 
'log mo 



mf(m I mo)(m / mo)^d log m 



r 

Ji( 



•log(mi+m2) 
'log mo 



f(m I mo)(m / mo)^d log m 



(16) 



The width of the Gaussian is chosen to obey two constraints: 

1. the Gaussian should have a sharp cut-off' at masses larger 
than the total mass of colliding aggregates. 

2. the power-law component must dominate the low mass part 
of the distribution. 

To obtain a representative fragment distribution, we average 
the mass spectrum over the impact parameters (see Sect. 14. lb . 
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Fig. 14: The filling factor averaged over the impact parameter vs the 
collision energy scaled with the number of monomers N and in units 
of the rolling energy E^^w. The dashed areas correspond to the spread 
around the mean. The dotted lines indicate the initial value of the filling 
factor for the colliding particles. 



Some impact-parameter-averaged distributions are shown in 
Fig. [131 The panels show the evolution of mass spectrum for 
both compact and ffuff'y aggregates. The eff'ect of impact averag- 
ing is immediately visible through the presence of the large frag- 
ment component at all energies. At low energies, growth domi- 
nates, and from central collisions one would expect to only see 
the large mass component at the total mass. However, grazing 
collisions contribute in the original, pre-impact mass bin, and 
reduce the average mass of the large fragment component. An 
increase in the energy leads to the onset of erosion. Small frag- 
ments appear with a very ffat slope of the power-law distribu- 
tion. Further increase in the energy enhances the erosion. The 
component of small fragments grows and begins to steepen the 
powerlaw slope slightly. Eventually, growth changes into frag- 
mentation. The large fragment component, however, remains in 
the distribution. Using impact-parameter-averages reduces both 
the eff'ect of fragmentation and diff'erences between ffuff'y and 
compact particles over what would be expected from head-on 
collisions. 



4.4. Compactness evolution 

The filling factor 0cr, similarly to mass spectrum, shows a diff'er- 
ent behavior in central and off'set collisions, respectively. The 
eff'ect of impact parameter averaging is presented in Fig.[T4l 
Compact aggregates are decompressed. At higher energies of 
about E = NE^oih aggregates undergo compression that mostly 
compensates the off'set collisions and the dimer-like structure 
discussed in Sect. 13.21 

Fluff'y aggregates also increase porosity in low energy colli- 
sions. The increasing energy, however, causes both a decrease in 
the mass of the largest fragment and compaction of the aggre- 
gates. Grazing collisions at a high energy produce aggregates 
with only weakly changed structure. Therefore, the structural 
change is dominated by compressing head-on impacts. A further 
increase in the collision energy results in stronger compression. 
The maximum filling factor of about (per ~ 0.19 is reached at 
E = NE^ow. Any further increase in the energy does not aff'ect 
the porosity. Particles are completely disrupted and the filling 
factor is dominated by small, ffuff'y fragments. 
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Fig. 13: Distributions averaged over the impact parameter. Results for fluffy and compact aggregates at energies A\ E - 5.7 • 10 ^NE\^^, B: 
E = 5J • lO-^TV^br, C:E = O.lSA^^br, D: £ = 0.927V£br, E: ^ = ^2NE^,, and F: £ = 14.6A^^br 



Small fragments produced by erosion or fragmentation are 
easily described by a single power-law. Regardless of initial 
porosity of the impact energy, small particles have the filling fac- 
tor given by Eq. ([T4l) . 

4.5. Format 

Our recipe provides parameters required to reconstruct mass dis- 
tributions averaged over impact parameter. The distribution of 
fragment masses is given by F{m) dlog m = f(m) dlog m, 
where f(m) dm provides the number of particles of mass m in 
mass interval between m and m -\- dm. Thus, the functional form 
of our recipe consists of the two components and is given by 

F(m) = ^im^ -\- —— — exp — , (17) 

y2ncrQ 2.(7 ^ 

where and ^2 are the normalization constants of the power- 
law component and the Gaussian, respectively. These constants 
are not provided since the distribution should be re-normalized 
in order to conserve the mass. Instead we give the mass ratio 
of the power-law and Gaussian components, which should be 
used to determine the mass and normalization constant in each 
component. The power-law component extends from a monomer 
mass to a quarter of the total mass. 

Our recipe is provided in tabulated form. The parameters 
required to reproduce the collisional outcome may be interpo- 
lated linearly. Our parameter space is covered very well and 
spans from very fluff'y aggregates of fractal dimension Z)f = 1.5 
through fluff'y fractal (Z)f = 2.0) and non fractal PC A aggregates 
to very compact particles of 0o- = 0.251. The energy space is 



also well sampled. Our recipe is based on simulations from hit- 
and- stick regime up to a catastrophic destruction. 

The main diff'erence between the local and the global sub- 
recipes stems from the mass ratio between impactor and tar- 
get. However, small projectiles, when carrying sufficient energy 
(i.e., impacting with very high speeds), may also shatter an entire 
large target aggregate, i.e. causing a global effect. Therefore, the 
global recipe must be used not only at mass ratio close to unity, 
but also at impact energies sufficient to globally affect the entire 
aggregate. Eigure [15] presents an algorithm used to distinguish 
between the local and the global recipes. When a collision of 
aggregates characterized by their masses m/, filling factors 0/, 
and some material properties, occurs at a given impact energy, 
one must check whether the collision is in the recipe domain 
(restructuring or fragmentation) or in the hit-and-stick regime. 
The test checks whether the energy is sufficient to cause any re- 
structuring. If the condition £^coii < 5£'roii is satisfied, the hit-and - 
stick mecha n ism is applied (for details see lOrmel et al.l (l2007l) . 
lOrmel et aP (l2009h ). Otherwise, a second check is performed. 
The global recipe is applied if any of the two following condi- 
tions is true: 

^>0.1, (18) 

m2 

or 

E > 5A^£br. (19) 

If none of the above is true, the local recipe must be used. 

The recipe provides a description of the outcome of a col- 
lision between aggregates of the same filling factor. To de- 
scribe collision of aggregates with different compactness, a mass 
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Fig. 15: Algorithm to choose between local and global recipes, 
weighted average filling factor must be determined 



(20) 



Therefore, the filling factor that should be used is dominated by 
that of the more massive aggregate. In particular, the outcome 
of a collision between particles of diff'erent masses depends very 
much on the porosity of the target aggregates and how deep it can 
be penetrated. Collisions of equal mass particles, on the other 
hand, are not dominated by one species. Therefore, the contribu- 
tion from both aggregates should be about equal. The compres- 
sion or decompression is expected to be weaker than for fluff'y 
particles, as compact aggregates are more resistant to restructur- 
ing (see Sect. [3]). 

4.6. A complete quantitative description 

The mass distribution of the collisional outcome for equal mass 
colliding projectiles can be constructed using parameters read 
from Tab. |7]-Tab. [6l Fig. fT6l shows contour plots of the required 
4 parameters. Intuitively, the mean mass of the large component 
(upper left panel) decreases with increasing energy. Similarly, 
the width of this component (upper right panel) decreases with 
increasing energy as a result of fragmentation and grazing colli- 
sions. Therefore, faster impacts cause formation of the largest 
fragments with the lower mass, and the tail of the Gaussian 
component decreases. This shattered mass shifts then to the 
small fragments power-law component. The power-law compo- 
nent does not exist at low energies and only large aggregates 




0.14 0.16 0.18 0.20 0.22 0,24 



Fig. 17: The change in the filling factor relative to the initial filling fac- 
tor cpf. cpjcpf 
compaction. 



are produced in this stage due to growth and grazing collisions. 
Thus, the slope of this component (lower left panel) starts to 
decrease only at larger energies. Note that for a weak erosion, 
where the slope of the power-law cannot be determined we as- 
sume without loss of generality q = 0. 

Fig- [13 shows compression of aggregates in the global 
recipe. This contour plot illustrates Tab.[9l The low-energy col- 
lisions cause a decrease in the filling factor 0^- for all aggregates. 
The eff'ect of decompaction is similar for all particles, although 
the strongest decrease in 0^- is observed for the most compact 
aggregates. Faster impacts result in compression and an increase 
in the filling factor. Aggregates with the lowest filling factor (pg- 
show the lowest structural change. The density increases only by 
a factor of about (po-/(pa^ = 1.1 1 at the highest energies. For very 
compact aggregates, at the same energy, this change is some- 
what higher (0o-/0J?^ = 1.22). This value is strongly aff'ected 
by fragmentation. The average largest collision remnant in this 
case consists of only 0.16 of the total mass. The maximum com- 
pression observed for head-on collisions at an energy of about 
E = N Eyow reaches only about 0.89 of the initial filling factor. 
This means that the average largest fragment is decompressed. 
Fluff'y aggregates, both fractal and non fractal, show similar be- 
havior. For these fluff'y particles the boundary between decom- 
paction and compression is at the energies of about E = N EyoW- 
The filling factor increases further at higher energies. The max- 
imum compression is reached at highest energies. The largest 
fragment at this energy is about a quarter of the total mass, ex- 
cept of the most fluff'y aggregates, that are very difficult to disrupt 
completely. 

4.6.1. Global recipe 

In this section we present tables containing the global recipe. 
They describe both the distribution of fragments (Tab. [S] Tab. [6l 
Tab. [71 and Tab. [8]) and the change in the filling factor of the 
largest fragment, relative to the initial (po- (Tab. [9]). 

The shape of the distribution of fragment masses depends on 
the number of contacts that can be broken. The higher the en- 
ergy is, the stronger and steeper the power-law component is. 
Therefore, the collision energy is normalized to the breaking en- 
ergy £^br. Moreover, the energy is distributed globally in the ag- 
gregates (hence the name global recipe). Therefore, the energy 
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Fig. 16: Parameters as a function of the filling factor and impact energy. Clockwise: mean mass of the Gaussian component, width of the Gaussian 
component, mass ratio of the power-law to Gaussian component, slope of the power-law. 



is also scaled with the total number of particles (this is ap- 
proximately equal to the number of contacts). The first column 
in each of tables [5] -[8] shows the normalized energy E/(N £^br). 
The following four columns indicate the recipe output quantity 
(see Tab.lU) for aggregates with diff'erent initial filling factors 
0JJ?\ Table [5]presents the slope q of the power-law component of 
small fragments. The mass normalization can be done with the 
total mass of colliding aggregates mi +m2 and the ratio Mr of the 
mass in this power-law component over the Gaussian component 
given in Tab. [S] 



Table 5: The slope of the power-law component. 



^/(A^^br) 


0.122 




0.155 


ni 

^0.189 


0.251 


5.72 . 10-4 


0.00 


0.00 


0.00 


0.00 


2.06 • 10-2 


0.00 


0.00 


0.00 


0.00 


5.72 • 10-2 


0.00 


0.00 


0.00 


0.00 


1.29-10-1 


0.00 


0.11 


0.00 


0.00 


2.29 • 10-1 


0.10 


0.12 


0.00 


0.00 


9.15-10-1 


-0.20 


-0.28 


-0.31 


-0.52 


3.66 • 10^ 


-0.69 


-0.49 


-0.42 


-0.53 


8.24 . 10^ 


-0.73 


-0.81 


-0.77 


-0.65 


1.46-101 


-0.85 


-0.89 


-1.03 


-0.93 



The Gaussian component can be reconstructed using three 
quantities; the ratio of the mass in the two components Mr, the 
mean mass Mg, and the width ctq. Table [7] provides the mean 
mass, and the width is given in Tab.[8l 



The change in porosity of aggregates depends on the number 
of contacts that can roll, as the rolling is the main mechanism 
responsible for the restructuring. Therefore, the energy in the 
first column of Tab. [9] is normalized to the rolling energy £'roii. 
The restructuring also aff'ects the structure of aggregates glob- 
ally. Thus, the collision energy is also normalized to the total 
number of monomers, which approximates the initial number of 
contacts. Table [9] shows the change in the geometrical filling fac- 
tor (per relative to the initial density 0cr^^^ This new filling factor 
represents large fragments from the Gaussian component. The 
power-law component is simply described by Eq. ([T4l) . 

The diff'erence in energy scaling applied to the part of the 
recipe describing the fragment distribution and to the part of 
the recipe providing the compaction, respectively, accounts for 
cases of diff'erent material properties, where the ratio of the 



Table 6: Mass ratio of the small fragments component (power-law) to 
the largest fragment component (the Gaussian) 





0.122 




0.155 


ni 

0.189 


0.251 


5.72 • 10-4 


0.00 


10- 


-u 


0.00 


10- 


-0 


0.00 


10- 


-0 


0.00 


10- 


-0 


2.06 • 10-2 


0.00 


10" 


-0 


5.43 


10- 


-5 


0.00 


10- 


-0 


0.00 


10- 


-0 


5.72 . 10-2 


1.17 


10- 


-3 


3.26 


10- 


-4 


8.14 


10- 


-5 


1.70 


10- 


-4 


1.29-10-1 


4.61 


10- 


-3 


1.47 


10- 


-2 


4.53 


10- 


-3 


2.86 


10- 


-3 


2.29 • 10-1 


3.67 


10- 


-2 


5.66 


10- 


-2 


1.14 


10- 


-2 


1.04 


10- 


-2 


9.15-10-1 


1.73 


10- 


-1 


2.67 


10- 


-1 


1.49 


10- 


-1 


1.39 


10- 


-1 


3.66 • 10^ 


1.91 


10- 


-1 


6.83 


10- 


-1 


6.63 


10- 


-1 


7.87 


10- 


-1 


8.24 • 10^ 


2.03 


10- 


-1 


7.53 


10- 


-1 


7.57 


10- 


-1 


1.72 


10- 


-0 


1.46- lOi 


2.02 


10- 


-1 


7.45 


10- 


-1 


8.79 


10- 


-1 


2.39 


10- 


-0 
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rolling energy £"1-011 to the breaking energy £br may be different 
than for Quartz. Therefore, different physical processes (com- 
paction/decompaction and erosion/fragmentation) are selected 
by different scalings, i.e., restructuring scales with the rolling en- 
ergy, while the fragmentation scales with the breaking energy. 

4.6.2. Local recipe 

The local recipe describes the outcome of a collision between 
a small aggregate and a large target. For high energy impacts, 
erosion occurs, resulting in what is basically the large target 
aggregate accompanied by a distribution of small fragments. 
The energy is locally distributed over monomers of the small 
impactor and some surface grains of the target. Thus, in this 
case we scale the energy with the reduced number of monomers 



Table 7: Mean mass of the Gaussian component normalized to the total 
mass. 



E/(NE^r) 


0.122 




0.155 


ni 

^0.189 


0.251 


5.72 • 10-^ 


0.671 


0.756 


0.749 


0.883 


2.06 • 10-2 


0.671 


0.756 


0.711 


0.883 


5.72 • 10-2 


0.671 


0.736 


0.711 


0.883 


1.29-10-1 


0.610 


0.697 


0.691 


0.883 


2.29 • 10-1 


0.575 


0.620 


0.671 


0.829 


9.15-10-1 


0.430 


0.436 


0.485 


0.536 


3.66 • 10^ 


0.427 


0.426 


0.423 


0.415 


8.24 • 10^ 


0.419 


0.418 


0.421 


0.401 


1.46-101 


0.422 


0.428 


0.436 


0.359 



Table 8: Width of the large fragments component (Gaussian) normal- 
ized to the total mass. 





0.122 




0.155 


ni 

^0.189 


0.251 


5.72 - 10-4 


0.124 


0.170 


0.118 


0.098 


2.06 - 10-2 


0.124 


0.170 


0.123 


0.098 


5.72 - 10-2 


0.158 


0.175 


0.123 


0.117 


1.29-10-1 


0.120 


0.169 


0.147 


0.137 


2.29 - 10-1 


0.141 


0.164 


0.170 


0.169 


9.15-10-1 


0.077 


0.107 


0.135 


0.164 


3.66 - 10^ 


0.061 


0.065 


0.061 


0.067 


8.24 - 10^ 


0.070 


0.066 


0.063 


0.071 


1.46-101 


0.068 


0.062 


0.061 


0.089 



Table 9: A fractional change in the geometrical filling factor (/)(j/(f)^^\ 
averaged over the impact parameter. 



EKNEroll) 


0.122 




0.155 


ni 

^0.189 


0.251 


8.44 - 10-4 


0.867 


0.818 


0.837 


0.796 


3.04 - 10-2 


0.874 


0.843 


0.860 


0.816 


8.44 - 10-2 


0.877 


0.867 


0.878 


0.822 


1.90-10-1 


0.923 


0.901 


0.902 


0.826 


3.37-10-1 


0.942 


0.944 


0.943 


0.838 


1.35- 10^ 


1.067 


1.053 


1.005 


0.888 


5.40 - 10^ 


1.082 


1.206 


1.144 


0.937 


1.21 - lOi 


1.088 


1.322 


1.274 


1.032 


2.16- lOi 


1.109 


1.368 


1.395 


1.223 



A/^ = N1N2KN1 +N2), which basically is the number of grains in 
the small aggregate. As the erosion is determined by the number 
of contacts that can be broken, we scale the collision energy to 
the breaking energy £br. 



Table 10: Mean ejected mass relative to the projectile mass. The mass 
ratio of the colliding particles is 0.001. 



E/(N^E^,) 










0.07 


0.09 


0.13 


0.16 


0.23 


0.97 


0.44 


0.52 


0.11 


0.92 


0.98 


0.58 


0.52 


0.11 


3.66 


1.15 


97.6 


3.42 


0.74 


14.6 


2.70 


109.3 


21.3 


7.83 


33.0 


1.21 


125.1 


30.9 


14.7 


58.6 


2.76 


148.0 


34.8 


26.9 



Table [To] shows the mean mass that is ejected during a col- 
lision, relative to the mass of a smaller aggregate. The first col- 
umn indicates the scaled energy, and the following four columns 
show the mean ejected mass for collisions of aggregates of dif- 
ferent initial filling factors (po-^^\ The mass of the cratered target 
aggregate can then be immediately calculated. 

The structure modification in the local recipe applies to the 
large target aggregate only, as the filling factor of small frag- 
ments follows a simple relation (see Eq. ([T4l)). Bombarding large 
aggregates with small projectiles results in very small relative 
change in the filling factor (that of the target particle), as the fill- 
ing factor quantifies the global structure of the large target aggre- 
gate. In this case, the energy scaling should thus be done in re- 
spect to the total number of monomers N, which is almost equal 
to the number of particles in the larger aggregate. Moreover, 
the restructuring mechanisms is determined by the number of 
monomers that can roll. Thus the energy is scaled by the rolling 
energy £:roii. 

Table [TT] presents the relative change in the geometrical fill- 
ing factor for the target particle. The energy listed in the first 
column is normalized to the rolling energy E^qw and to the total 
number of monomers N. 



Table 11: Relative change in the filling factor (pcr/(p^^^- Mass ratio of the 
two colliding particle is 10~^. 



E/(NE,oii) 


0.07 


0; 

0.09 


ni 

0.13 


0.16 


0.00034 


1.00004 


1.00001 


1.00005 


1.00031 


0.00135 


1.00002 


1.00011 


1.00964 


1.00072 


0.00539 


1.00010 


1.04656 


1.00983 


1.00051 


0.02155 


1.00077 


1.05293 


1.01031 


0.99677 


0.04849 


1.00012 


1.06677 


1.01024 


0.99314 


0.08620 


1.00082 


1.06368 


1.00689 


0.98781 



5. Conclusions and future work 

In this work we present results of the extensive parameter study 
of collisions of three-dimensional aggregates. The outcome of 
a collision is provided in terms of the mass distribution of 
fragments as well as the structure of the produced particles. 
These simulations agree with the experimental results in this 
size regime and provide scaling that allows for extrapolation to 
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slightly larger sizes as well as different material properties of 
monomers (i.e., composition and size of individual grains). 

Our simulations indicate new important mechanisms that in- 
fluence both the structure and the mass distribution of aggre- 
gates: 

- The restructuring of aggregates depends mainly on the col- 
lision energy. The compaction is reserved for head-on im- 
pacts, while off'set collisions produce elongated and decom- 
pacted particles. 

- In the case of erosion and shattering, the structure of 
small fragments can be written in a very simple form (see 
Eq. ([T4l)). regardless of the initial compactness of colliding 
aggregates or the impact energy. 

- The mass distribution of particles produced in a collision 
consists of two individual components. The power-law dis- 
tribution of small fragments is accompanied by a strongly 
pronounced component of large fragments. 

- The shape of the mass distribution for the collision outcome 
averaged over the impact parameter is generally independent 
of the impact energy for slow collisions. At high energies, 
however, the shape of the distribution is almost independent 
of the structure of the colliding particles. 

These points are combined together in a form of the recipe that 
provides the quantitative outcome of a collision. This recipe is 
formulated in a simple way that can be easily applied to mod- 
els of dust coagulation in various environments (e.g., molecular 
clouds or protoplanetary disks). In fact, this recipe has already 
been applie d to study growth o f small dust aggregates in molec- 
ular clouds (lOrmel et al.ll2Q09h . 
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